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' i A complete description of W and Z boson production at high-energy hadron 

Qj I colliders requires the resummation of large Sudakov double logarithms which 

r^ • dominate the transverse momentum (qx) distribution at small qx- We compare 

different prescriptions for performing this resummation, in particular implicit 

impact parameter space resummation versus explicit transverse momentum space 
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^ • resummation. We argue that the latter method can be formulated so as to retain 



the advantages of the former, while at the same time allowing a smooth transition 
to finite order dominance at high qx- 
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1 Introduction 

The production of W and Z bosons at hadron collliders provides several fundamental 
tests of perturbative QCD. A problem of particular topical interest is the calculation of 
the transverse momentum {qt) distribution of the produced vector boson. Data from 
the Tevatron pp collider experiments now cover the regions of both small and large qx 
||l]|. For large g^, fixed-order perturbative calculations should be sufficient, and indeed 
the current 0(a|) predictions agree well with experiment. However, it is the small 
qt region which is more theoretically challenging as here one encounters the infrared 
structure of the theory in the form of large logarithms of Q/qt, where Q is of the 
order of the weak boson mass . The presence of higher-order contributions of order 
agin "" (Q'^/Qt) leads to a breakdown of fixed-order perturbation theory as qt -^ 0. 
Although there is much data available, a completely consistent theoretical treatment 
in agreement with experiment has not yet been developed.^ Furthermore the most so- 
phisticated treatments ||^, involving the resummation of the large logarithms in impact 
parameter space, lead to results that are in practice difficult to merge with the large qr 
fixed-order expressions. The form of the low qt distribution is not only of theoretical 
interest - for example, a proper description is needed for an accurate determination 
of the W boson mass. The formalism also applies directly to the production of any 
massive colour-neutral particle, including the Higgs boson. 

The purpose of this paper is to investigate the resummation of large ln(-^) log- 
arithms directly in momentum space. The question has been addressed recently in 
Refs. 0] and 0, and our analysis can be considered an extension of these studies. 

We begin with a simple formulation of the problem. The large logarithms in the W 
transverse momentum distribution at small Qt arise from soft gluon emission from the 
incoming (massless) quarks in the basic qq' — >■ W process. The logarithmic dependence 
enters the differential cross section formula through terms proportional to (for J^) 

agln'"-"^(^] m = l,...,2n. (1) 

In practice, the coefficients of the terms represented by the expressions in (|^) are known 
only for m = 1,2, 3, 4; the rest remain unknown. In the simplest case, the cross section 
at small qx is approximated by summing only the leading logarithm terms in (|I|) (i.e. 
m = 1), which gives a Sudakov form factor (for the qq' -^ W subprocess cross section): 
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^For a review of the literature, see for example El. 

^For simplicity we will only refer to the W boson in our study, but obviously all our conclusions 
apply equally well to Z production also. 



This double leading logarithm approximation (DLLA) corresponds to the situation 
when only the contributions of soft and collinear gluons are included and the strong 
ordering of the gluons' momenta is additionally imposed 0. The fact that the Sudakov 
form factor in (0) results in a suppression of the cross section as g^ — >■ indicates that 
sub-leading logarithms (i.e. m > 2 in (|1|)) have to be taken into account in the 
resummation. Two methods of doing this have been proposed. 

(i) Resummation in impact parameter h space (where 6 is a two-dimensional Fourier 
conjugate of gr)- This method correctly takes into account the known logarithmic 
terms and also certain kinematic features of the gluon emission. However it 
suffers from several drawbacks, see [|, ^, in particular it leads to an unphysical 
extrapolation to large g^, as will be illustrated below. 

(ii) Resummation in g^-space (proposed recently in Ref. 0). This method in- 
cludes exactly the same logarithms as the 6-space formalism up to and includ- 
ing agin "~ (%) terms, but omits sub-leading 'kinematic' logarithms, i.e. sub- 
leading logarithms which arise when the strong ordering assumption is relaxed 
and which are automatically included in fe-space. These terms start to contribute 
from the fourth 'tower' of logarithms down onwards. The question is whether 
it is possible to include sufficient kinematic logarithms using this technique that 
the 6-space cross section can be adequately approximated by resummation in g^ 
space in the region of qt relevant to the comparison with data. Furthermore, in 
this approach the problems with merging the resummed and fixed-order results 
in the 'medium' qt region can be more easily circumvented. 

In this paper we explore further the g^-space resummation approach. We extend the 
work of by including all NNNL logarithms and higher 'towers' of known kinematic 
logarithms. The effect on the cross section of systematically adding sub-leading terms 
is quantified. The goal is to achieve a momentum-space-resummed cross section which 
reproduces the impact-parameter-resummed cross section in the regions of qt relevant 
to the experimental measurements and which includes all known calculated coefficients. 
The paper is organized as follows. In the following section we present the basic 
theoretical expressions for the resummed cross section in both b and qx spaces. In 
Section |] we consider in turn the quantitative effect of the various sub-leading logarithm 
contributions, in particular those from higher-order coefficients, kinematics, and the 
running of as- Section ^ contains a summary and some conclusions. 

2 Theoretical description of the small Qt distribu- 
tion: resummation in qt space 

The general 6-space expression for the differential cross section for vector boson pro- 
duction at small qt has been given in M (and see references therein). For simplicity. 



throughout this study we shall restrict our attention to the parton-level subprocess 
cross section: parton distribution functions can in principle be incorporated to yield 
the hadron level cross section without changing any of our conclusions on resumma- 
tion.0 Using the 6-space formalism, at the purely perturbative level, one obtains for 
the differential cross section [H| 
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and Q = Mw in the present context. 

The first two coefficients in each of the series in (|^), i.e. A^^\ A^'^\ B^^\ B^'^\ 
are well known P|. They can be obtained from the exact LO + NLO perturbative 
calculations in the high qt region by comparing the logarithmic terms therein with 
the corresponding logarithms generated by the first three terms of the expansion of 
exp(S(6,g2)) in d). Explicitly, 
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with Cp = 4/3, Tr = 1/2 and A^ = 3. It is instructive to see how the logarithms 
in 6-space generate logarithms in g^—space. For illustration, we take only the leading 
coefficient A^^^ = 2Cf to be non-zero in e'^^'*'*^ \ and assume a fixed coupling as- This 
corresponds to 

da an f^ ^ „ ^ , ,, / a.^Cp , fQ'^h'^^ 
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■^Formally, one can imagine taking (t — Q'^/s)^ moments of the hadron cross section, which aUows 
the subprocess cross section to be factored out. Our subprocess cross section corresponds to the iV = 
moment in this sense, ^ = ^S(0, ^, ^,as{n'^)), see Ref. §. 
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The expressions are made more compact by defining new variables i] = %, z = b^Q"^, 
X = s^^ zo = 4exp(-27s) = bl Then 

(To drj 4 Jo 

and we encounter the same expression as in 0] , which describes the emission of soft and 
colhnear gluons with transverse momentum conservation taken into account. The resuh 
of numerically integrating (|^) and its comparison with the DLLA approximation (^ 
is shown in Fig. |I|.0 The cross sections are similar over a broad range in t]: the main 
differences occur at (i) small rj, where the DLLA curve is suppressed to zero and the 
6-space curve tends to a finite value, and (ii) at large rj (strictly, outside the domain of 
validity of either expression). In the latter case, the g^-space cross section vanishes at 
?7 = 1 by virtue of the overall factor of ln(l/?7). This is a crude approximation to the 
(formally correct) vanishing of the leading-order cross section at the kinematic limit 
~ Mw- In contrast, the 6-space cross section has no information about this kinematic 
limit, and is non-zero at r/ = 1. Furthermore, at large qt the 6-space cross section 
oscillates about 0. This can be seen in Fig. |^, which extends the cross section of Fig. [^ 
to large rj on a. linear scale. The first zero of the oscillation is clearly evident. Now since 
this occurs far outside the physical region it might be argued that it is not a problem 
in practice. However, when the first sub-leading logarithm B^^^ is included, the first 
zero moves inside the physical region, as shown in the figure. It is this behaviour 
which causes problems in merging the large qt fixed-order result with the resummed 
expression, since the compensating terms have also to be given an unphysical oscillating 
form. 

After performing a partial integration and expanding terms under the integral in 
(|^), one obtains 

which is now a perturbation series in q^ space. Here the numbers 6^(00) are defined 
by 

b^{oo)= rdyMy)W'^{^) (9) 
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and can be calculated explicitly from the generating function 
E;^«„.(oo)^exp[-2EL)(-) 

so that e.g. bo{oo) = 1, bi{oo) = 62(00) = 0, 63(00) = —^({3) etc. 

^For all our studies with a fixed value of the coupling we take A = "^ -^ = 0.085, n,/ — 4. 



Before studying the various approximations to the cross section in detail, it is 
worthwhile comparing the two expressions for the differential cross section (|^) and (§). 
The first, in b space, calculates the cross section in terms of a one-dimensional integral. 
The cross section is well defined at all values of qt, and in particular at qt = 0. In 
practice, however, non-perturbative (e.g. intrinsic /cr smearing) effects dominate in this 
region. These can be modelled by introducing an additional large-6 suppressing piece 
in the integrand, for example exp(— a^fe^) where a is a measure of the kx smearing. In 
practice, since the argument of the running coupling in 6-space is proportional to 1/6, 
some form of large-6 cut-off or freezing must also be performed. 

As already mentioned, however, the large qt behaviour of (|^) is not physical: 
the integral has no knowledge of the exact kinematic upper limit on qt, although 
numerically it becomes small when qt ~ Q- More problematically, as qt is increased 
the distribution starts to oscillate, and it is this feature (built-in via the Bessel function) 
which makes it difficult to merge with the finite-order large qt cross section. 

In contrast, the qT-spa.ce cross section (^ is an asymptotic series. The logarithms 
are singular at gr = 0, although as argued above this is in any case the region where 
non-perturbative effects dominate. As with any asymptotic series, care must be taken 
with the number of terms retained. Merging with the fixed-order large qx cross section 
is in principle straightforward, since the g^ logarithms of (|^) can simply be removed 
from the finite-order pieces to avoid double counting. 

3 Quantitative study of resummed cross sections 

The sub- leading corrections to (0) have three origins. First, there are sub- leading terms 
arising from the matrix elements that are associated with the coefficients A^'^\ B^^\ etc. 
Secondly, there are sub-leading effects resulting from the running of the strong coupling 
as- Finally, there are also sub-leading terms in the form of kinematic logarithms, 
always appearing with 6m (oo) (m > 1) coefficients. In the following subsection we 
focus on this particular type of sub-leading effect and assess its importance. We isolate 
the effects induced by kinematic logarithms by fixing the coupling and taking only 
leading terms arising from the matrix element. Then, in the following subsections, we 
subsequently switch on running coupling effects and sub-leading logarithms from the 
matrix element. 

3.1 Fixed coupling analysis 

We begin our study of (H) by performing the resummation for the simple case m = 0. 
This is the DLLA of Eq. (Q), i.e. all radiated gluons are soft and collinear with 
strong transverse momentum and energy ordering, and no account taken of transverse 



momentuni conservation: 

i?-iEA"tiZH)^M=o).n-(i). (U) 

aodr] 7]^^ (A^-lj! \r]J 

Next we investigate the effect of including all the m > terms in @. (To be 
approximated numerically the series has to be truncated at some iVmax- Thus full eval- 
uation of (|]) up to the A^max-th term requires knowledge of the first 2A^max— 1 coefficients 
6m (oo).) The first 20 coefficients, calculated according to (|10|), are listed in Table |1|. We 



find (see Fig. ^ that for large m the coefficients behave as b„i{oo) ^ C {—l)^m\2~ 
where C is a constant. 

Taking more terms into account, i.e. m > 1, we find that, as expected, the sum (^ 
exhibits behaviour consistent with an asymptotic series. A single (A^, m) contribution 
to (H) is of the form 

and to show the complexity of the resummation (^) we display these individual con- 
tributions in Fig. |[ For all rj the biggest contributions arise when m ~ 2A^ — 1, since 
the coefficients b„i{oo) are largest there. As t] decreases, contributions with smaller m 

2N~m~l (\ 



become more significant due to the terms In 

Our first task is to investigate numerically the dependence of (^ on the point of 
truncation A^max, i-e. the order of the perturbative expansion, and the number of terms 
included in the internal summation (^ - the 'cut-off' value mmax, equivalent to the 
number of known 'towers' of logs down from leading. Obviously for different pairs 
(^max^w^max) different contributions (|T2p are summed, see Fig. |^. In Fig. ^ we show 
a 3D cumulative plot of (||) which illustrates some of the features discussed below. 
Each plotted value for a given point (A^max, 'Timax) represents the sum (H), truncated at 
A^max and calculated with 6m (oo) = for m > mmax- The distinctive plateau present 
for large values of A^max and small mmax is equivalent to recovering the 6-space result 
for various 77. Notice how for smaller ri this plateau has a tendency to contract. If all 
2A^max — 1 = ""^max Coefficients 6m (00) are taken into account (see Fig. |a), it seems that 
the 6-space result cannot be approximated for any value of A^max, except for the region 
of large t]. This should not be surprising, considering that the 'towers' of logarithms 
have been truncated. Conversely, if only the first few coefficients (mmax < 2A^max — 1) 
are known (Fig. ^o) and the rest of them are set to zero, then in some sense one is closer 
to the DLLA situation and it is possible to find A^max such that the cross section (H) 
approaches the 6-space result, at least for the range of r] considered here. Moreover, 
it can be seen from Fig. | that it is necessary to consider the first jew coefficients to 
achieve the best approximation of the 6-space result. 

So far we have not attempted any analytic resummation of the series for the 
gy-space cross section given in (^). It is interesting to see whether factorizing out 

6 



the resumined DLLA piece from (P) leads to an improvement in the approximation of 
the 6-space result. We again start from (|) but now we extract the Sudakov factor 
exp ( — I In^ (-)) from the sum to get 
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A key feature of ([T3|) is that after extracting the Sudakov factor, the residual perturba- 
tion series has at most N+ 1 logarithms of 1/t] at A^th order in perturbation theory, i.e. 
the leading terms are now A^ In "*" ( - j . However we know that these terms must sum 
to give a large contribution as r/ ^ in order to compensate the overall suppression 
from the Sudakov factor. 

The terms which contribute to the new series (0) are illustrated schematically in 
Fig. H. Notice that the extraction of the Sudakov factor results in an ability to sum 
an infinite subseries of logarithms. This observation constitutes a basis for our further 
analysis. However, there is a shortcoming: in order to sum the first m 'towers' fully 
we need to take A^max = fn which leads us to include extra sub-leading contributions 
from more than the first m 'towers', cf. Fig. Rib. 



The individual contributions to the summation (\[^) 
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are presented in Fig. 0. As was the case for (|T2|), the importance of the In (-) factors 
diminishes for small m and large rj. Note that for small rj the sizes of the contributions 
are much smaller here than for (^). 

Next we perform the same analysis as for (|]), i.e. we study the dependence of (|T3D 
on A^max and on the 'cut-off' value rrimax- Comparing the cumulative plot for ([l^), 
Fig. ^, with Fig. |], we note that the 6-space result is now better approached close 
to the line m-max = ^max — 1 (such an effect can be observed in the case of (|) only 
for large rj). When mmax < ^max — 1 the balance between different contributions is 
apparently spoiled until mmax becomes considerably small. Again, it turns out that it 
is necessary to incorporate the first few sub-leading kinematic logarithmic terms (i.e. 
some moderate A'max) to obtain the best approximation of the 6-space result. 

The asymptotic properties of (0) can be most easily seen for large values of t] 
and large A'max- The apparent discrepancy between the 6-space result and the sum- 
mation ( pi^D is caused here by contributions with m = N — 1, i.e. terms propor- 
tional to &2Af-i(oo). The other terms are negligible as they contain powers of In f -j 
which for these r] are very small. The dominant contribution is then proportional to 



(^-l)^'-^^^^-^2^'^''''b2Nrr,^^-i{oo) / {N^^a.^- 1)\, i.e. the sign varies as (-l)^'"^". Fortunately, 
the range of t] for which the discrepancy occurs is outside the region of interest for the 
present discussion. Also, in practice we never take A^'^ax so large as to make this effect 
substantial. Nevertheless, it emphasises the necessity of performing a careful matching 
with the fixed-order perturbative result. 

We may therefore conclude that the expression (pTSf), with the Sudakov factor re- 
summed and factored out, enables us to resum infinite series of logarithms while at the 
same time allows for the inclusion of kinematical sub-leading logarithms. Moreover, it 
has very good convergence properties over a large range of r], while summing leading 
and sub-leading logarithmic terms. It is thus well suited for the purposes of this anal- 
ysis, i.e. reproducing the 6-space result by explicit resummation in qx space, and we 
will continue to use it for the rest of this study. 

3.2 Running coupling analysis 

The ultimate goal of the work presented here is to obtain a more accurate description, 
if possible, of the W production distribution. To this end, one has to incorporate 
various other sub-leading effects in addition to the kinematic logarithms discussed in 
the previous subsection. One example is the incorporation of the running of the strong 
coupling as into the cross section expression. This is achieved by substituting 

+ ^^ H^2) U In(^) - a) I + 0(4(/^^)) , (15) 
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/5o = 4 (ll - 3^^/) , (3r = - [102 - -Ufj , (16) 

in (plsl) . The effect on the DLLA form factor is to introduce a sequence of sub- leading 
logarithms whose coefficients depend on the /5-function coefficients defined in (|T3]). If 
only the 1-loop running of as is introduced then the cross section expression (|I3D reads 
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Notice the appearance in (0) of sub-leading 0(a|ln^(l/?7)) terms in the exponent 
of the Sudakov form factor. Interestingly, these logarithms can be eliminated by a 
particular scale choice: /i^ = Q'^'rf''^. This restores the same form as in ([T3| ) i.e. 



A(i) 



exp ( — "^^" In^ (-)); but now with a coupling which also depends on t^.Q However 
not all /^Q-dependent logarithms are eliminated. For example, we can see from (|18D 
that corrections of order /3o«5 In (-) remain. Of course in a complete calculation the 
dependence on the scale choice should disappear. To illustrate the residual dependence 
on the scale of the cross section ([T7|) we show (Fig. 10) results for two different choices: 
/i^ = T^Q^ = g|. and /x^ = Q'^rf'^. Also shown is the effect of truncating the sum in 
(P^ ) at first, second and third order .0 

We see that with ji^ = Q"^ rf/^ there is slightly less stability with respect to the 
truncation point than there is with the choice /i^ = q^. Note also that now the 



plots of ^1^ are no longer explicitly independent of Q, as was the case when as was 
taken to be constant. The additional dependence on Q enters into (0) via the coupling 

o^silJ^) = ct-siQ'^^^Qr )• Fo'^ values of qt where perturbative QCD can safely be applied 
(e.g. gT > 3 GeV) and the energies considered here, the resulting scale /i is bigger than 
the b quark mass, and we therefore avoid evolving the coupling through any quark mass 
thresholds. 

An extension of this analysis for the case of the two-loop running coupling is 
straightforward and gives 



J_da 
(TO dri 



as{f^')AW 




N-l 



N-1 



(A^-l)!^ 



N 



m 



m — k 
I 



I m k I N~m~k~ 
^2 1-31-41-1 



l-l 



X 



4:C4bN+m+2k+3l+2{00) + SCsfo 



'Ar+m+2fc+3/+l(oo) + 2C2&Af+m+2fc+3K°°) 



+Clb]\f+rn+2k+3l-l (oo) 



(19) 



Here 



On 



An 



InM- 



1 + ^ (/5o«5(/x') + ^AaK/u')) (2 In (i) - 3 In (^% 



^There is an analogous ycut~dependent scale choice for resummed jet cross sections in e+e^ anni- 
hilation, see for example H. 

^Obviously, now we truncate the sum in ( p^ ) at a certain order of asin'^), so that W in Fig. || 
should be read as 'power of as'- 
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Now, after fixing the choice of scale as described above, a new term proportional to 
I3la\{(j?) In^ (-) appears in the exponential. As a result, the line of points correspond- 
ing to terms of the form a^ In ~ (-) in Fig. ^ja will now change to the set of points 
corresponding to a|^ In''^ (-) or higher powers of In {-). Unlike the In^ (-) terms, those 

with In (-) do not cancel when /x^ = Q'^'^qr/ is chosen. It is impossible to cancel these 
terms by any choice of the renormalization scale. On the other hand, the choice which 
eliminates those terms with In (-) is also the one that minimizes the coefficient of the 
In*^ (-) terms. 



3.3 Resummation including all types of sub- leading logarithms 

A derivation of the expression for the cross section for the case when all known leading 
and sub-leading coefficients, i.e. y4*-^\ B^^\ A^'^\ B^^\ are taken into account follows the 
method introduced above and gives (fixed coupling to begin with) 
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Now each logarithmic term in the sum (^0]) acquires a factor which is a combination 
of the As, Bs and 6(00) 's. Notice that ahhough the various sub- leading logarithms 
are mixed together, they have a distinctive origin. We have mentioned already that 
the DLLA (i.e. retaining only terms of the form A^^^a^ln ~ (-) ) corresponds to 
the situation where all gluons are soft and collinear and where strong ordering of the 
transverse momenta and energies is imposed. We also know that other terms with A^^' 
multiplied by the bm{oo) arise from using the soft and collinear approximation for the 
matrix element but relaxing the strong-ordering condition. The sub-leading terms with 
B^^' etc. correspond to the situation where at least one gluon is either non-coUinear 
or energetic. 

In addition, extracting a Sudakov form factor from the sum (|20|) 'squeezes' it down 
to a summation over m from to A^ — 1, thereby reducing the number of fully known 
'towers' of logarithms: from the first four to the first two. This can be appreciated by 
comparing the logarithmic coefficients appearing in the expansion of the cross section 
up to and including the first three orders in as with (Table ^ and without (Table ^ 
the Sudakov form factor extracted. 

These tables also reveal another relevant property: sub-leading coefficients like 
A^"^^ are associated with b(oo) factors whose indices are not as high as those which 
accompany A^^\ Both In (-) terms and 6(00) factors originate from the same In (^2-) 
terms appearing at the very early stage of the derivation. In particular, let us focus 
on a term with a particular power of 0:5 and In (i). To reproduce such a term one has 

to take In (^ji-) up to the appropriate power, depending on its associated coefficient 

"0 

(i.e. A or B). For sub-leading coefficients this power will be obviously lower than for 
more leading ones. Hence the indices of the corresponding 6(00) factors are lower for 
more sub-leading coefficients. For example, the index of b{oo) accompanying the first 
sub-leading unknown coefficient A^^^ would be at least four less than the index of a 
corresponding 6(00) for the A^^^ coefficient, for given powers of as and ln(i). This 
observation provides us with a strong argument for justifying the inclusion of known 
parts of logarithmic terms from sub-leading 'towers' lower than NNNL. Physically, the 
b{oo) factors are of kinematic origin and as such they are much more relevant to the 
cross section than perturbative lower order coefficients in the expansion of Alasi^"^)) 
and B{as{iJ^))- While, as we have shown earlier, resummation of the terms containing 
6(00) factors with higher indices can still be of numerical importance for the final 
result, terms with unknown higher-order perturbative coefficients seem to contribute 
corrections of much smaller size. 
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A comparison of effects induced by the inclusion of successive sub-leading coef- 



ficients in (20) is demonstrated in Fig. 11. Also shown are the results of an exact 



integration in 6-space. The agreement for low values of r] is excellent. At high r] we 
encounter a discrepancy of the same nature as that discussed for the case of the lead- 
ing coefficient A^^\ i.e. for large A^max the expression (pOD either grows significantly or 
acquires negative values, depending on the number of terms at which the expression is 
truncated. For the values of A^max we use in our calculations this behaviour does not 
occur .[] 

With the above prescriptions for the sub-leading logarithmic terms, kinematic ef- 
fects and running coupling, we finally obtain a 'complete' expression for the cross 
section: 



J_da_ 
ao drj 



c^s{f^')A^'^ 5„^ /-ac,(/x2)AW 





Af-l 



IT 



A— m — 1\ Y^ i\ — m — k — 1 




N-m-k-l-l\ ^-™-^'-J-i / N _m-k-l-]-l 



/ i=0 

6 



m k I J i iV— m— fe— ( — J— I— 1 V^ h / \ /oi \ 

X C2 C^C^d^CQCi 2^ nCnbN+m+2k+3l+4:j+5i+n-2{00) . (21) 



n=l 



The expressions for 5*^ and the q are now of course much more complicated. Explicit 
expressions up to fifth order in the coupling constant 0:5 (/i^) are presented in the Ap- 
pendix. The fifth order appears here as a consequence of using the two-loop expansion 
of the running coupling and including the first four terms in the expansion of Sr^ (^. 
Numerical results based on the complete expression are displayed in Figs. IT^-ITB 



for the scale choice fi^ = Q'^^^qr/ . First, Fig. O shows the cross section with all 



four Ai, Bi coefficients included, together with the first three orders in as{^i'^) in the 
entire sum over A^ in (|21|). Notice the rapid convergence when the higher-orders are 



included. Fig. |13| shows how the first-order cross section is infiuenced by the various 
Ai, Bi coefficients. Notice how the relative impact on the leading-order Ai result 
is essentially Bi > A2 > B2, as might be expected. The effect of including B2 is 
numerically small, indicating a reasonable degree of convergence from the higher-order 
coefficients. 

The result (^1]) is valid in general for any choice of the renormalization scale fi. 
However, the expressions for S"^ and the c's will change depending on the particu- 
lar scale choice. For /i^ = Q^ the expression for the cross section -^^ corresponds 

to F^'^'^^qt) in Ref. [^. More precisely, F^'^'^^qx) as defined in |^ can be obtained 
^This is also true for the case of the running coupUng. 
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from (pTD by choosing an upper limit of summation A^ = 5 and putting 6j(oo) = 
for i > 1. This recipe comes from the observation that the only coefficient appearing 
together with bo{oo) is Ci, which contains terms up to C(a|((5^)).0 Moreover, the ex- 
pression (0), when calculated for /z^ = Q^ and expanded in powers of 05 (Q^) gives the 
0{a'g{Q'^)) 'perturbation theory' result F^P\qT) from ||^. In Fig. |T^ we show F^^^qx) 
and F^'''^\qT) as functions of gr, analogous to Fig. 3 in f^.^ 

The g-r-space formalism as presented in ||^ does not take account of non-zero bi{oo) 
with i > 1, and therefore does not include known contributions from sub- leading 
NNNL 'towers' of (kinematic) logarithms. This is partially compensated for in [^] by 
a redefinition of B^'^\ 

^(2) _ ^(2) ^ 2A(i)C(3) (22) 

which, although it does correctly account for a 0{a'g{Q'^)) term from the NNNL 'tower' 
(see Table ^, distorts other terms from this 'tower '.p^ With the help of our expres- 



sion (21) we can obtain the first four 'towers' fully resummed. It should however be 
remembered that the result is not 'pure', in the sense that it contains additional sub- 
leading terms. On the other hand, an effect induced by the redefinition of B^'^^ seems 
to be comparable with that caused by summing these sub-leading terms in (pi]) (with 
non-zero 63(00), 65(00)), see Fig. |15[ A difference obviously arises when the fourth 
tower is also resummed (cf. Fig. [l^) — numerically we encounter an increase in the 
cross section of approximately 3% for all values of qx, when the scale equals jj? = Q"^ . 
Furthermore, Fig. [l^ shows the effect of including 3,5,6,7,8 'towers', normalised to the 
4th-tower result, now for the scale choice jj? = Q'^^^q^ . The figure clearly shows the 
numerical importance of the kinematical logarithms of the higher towers, and also the 
stability over a broad range of relevant qx as more towers are included. However, this 
change is approximately of the same magnitude as the one observed for the change of 
the renormalization scale. In particular, for the 4 towers of logarithms, changing the 
scale from our default /x^ = Q'^^^qJ to the (lower) scale /x^ = q^ and the (higher) scale 
/i^ = Q^ changes the qr distribution by less than ±3% over the complete low-gT range. 

4 Summary and conclusions 

The g^-space formalism for describing vector boson production in hadron collisions is 
known to overcome many of the problems faced by the 6-space method. In this paper 
we have further investigated the gr-space approach. For the parton-level subprocess 



^Notc that throughout this work we consider the two-loop expansion of the running couphng. 

^Although we agree with [p| in the analytical result for the perturbative expansion, there is a 
significant numerical discrepancy that we are unable to account for. 

^° We also find a disagreement with the curves in Fig. 4 of Ref. [^. Since S'^' > i?^^-' > 0, after 
the replacement of B^^-* by B^^-* one would expect a smaller value for the cross section, in contrast to 
the displayed result in Ref. [g. 



13 



cross section we have modified the existing approach in order to incorporate sub- 
leading 'kinematic' logarithms. We have carefully studied the effect of various sub- 
leading contributions: higher-order perturbative coefficients, the running coupling, and 
the 'kinematic' logarithms. We have confirmed that the 'kinematic' logarithms are 
particularly important at small qt-, where they serve to cancel the suppressing effect 
of the Sudakov form factor. 

Our technique enables us to resum the first four logarithmic 'towers' including the 
NNNL series, together with the first few sub-leading 'kinematic' logarithms. We have 
shown that that the most significant quantitative change in the predictions for the 
cross section is caused by resumming the NNNL 'tower'. The fact that the fourth 
'tower' only changes the cross section by about 3% shows that the convergence of our 
expansion is certainly adequate for phenomenological applications. We note that a 
drawback of this method is an inability to select a particular numbers of 'towers' to be 
fully resummed. 

In this paper we have concentrated only on the perturbative contributions to the 
cross section. It is well known that in practice non-perturbative effects ('fcT smearing') 
are also important, and affect the W qt distribution in the very low qt region, see for 
example Ref. Q. These must be taken into account before assessing the impact of the 
sub-leading logarithmic contributions on the physical cross section. We will address 
these issues in a forthcoming study. 

Note added: As this paper was nearing completion we became aware of an interesting 
new study of soft gluon resummation, Ref. |]lOl, which addresses the same problem. 



In Ref. [0], resummation in 6-space is studied, in particular the impact of factorially 
growing terms in the a^ expansion therein. A closed expression is subsequently ob- 
tained for the corresponding qt space result, which resums logarithms at the NLL level 
in the Sudakov exponent. The kinematic logarithms are also treated differently, such 
that a singularity is encountered as qx — ^ q'^^^ > 0. We will compare the phenomeno- 
logical implications of the two approaches in a forthcoming study. 
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Appendix 

In this appendix we list the expressions for 5*^ and the q coefficients in (0), for the 
choice of the renormalization scale /i^ = Q'^r]^^^ . 



Sr, - --A 



1 «.(/.^)^/?o^A(^)ln(;ij , Po' asif^r Ai^^ 



270 vrMW V 135 7rM« 



V540 TT^AW 1620 7rM(i) / ^^^ ' J \r] 

f 1_ PlasifJ^YAi^ _ ( 1 AM(^) 5 /3o^Ai^(^) \ ,5 

^U2 vrMW \ 144 7rM«^216 vrMW j "^^^ ^ 

1 /3oAA(^) , 5 P,^B^^) \ ,, iplas{ilf\ (I 

"72 ttMW ^ 216 ttMW y "^^^ ^ ^ 18 tt^ r\r^ 



^ -1 (5,' asili^f B^^) ^ 2 P,P,as{li'YB^^) ^ 1 ^^ asji^-'f B^^) ^ 2 as{li'? P,' B^') 



9 vrMW 9 vrMW 3 tt^ AW 9 7r2 AW 

15 



X ln( - 



1 P^asif^^lB^^ ^ fiPoB^^ip^B^\ ,2,2 



vr 



-An 



1 A(2) 1 /?ofi«\ 



12 7r3A(i) 12 7r3A(i), 



ttMW 12 7rM( 



^)«,(^^) + l)ln(-) +( 



TT 



^(1) 



^I«J'" 



Cl 



, ff^B'^^nhy , (f&^^lSL^yM 



TT 



5^(1) 



30 



^.j 



^ 30 



1 (^27-4n^ + T AW )^^[v) 1 ln(i)5/5oM(2) 



y]-U 



162 n^AW 



asifi 



2\5 



/ 



+ 



^ 30 



^10 7r/303B(2) I 10 nj3oJhAW_\ ,^ /^l^ \ 



T AW I^^T^ 



9 ttMW 



27 vr^yK^) 



TT"' 



V 



«s(At 



2\4 



+ 



lA5(^)ln(i) ^ (s^^l^^ + lO^^^^Vn^i 



3 7r3A(i) 



^ 30 



A(i) 



n- 



6 7r3 A(i) 



i^^\ 



V 



^^(/i 



2n3 



+ is 



Ad) 



10 "'?n^'"' +10 ""^.Vr' lln(i 



A(i) j^'^U 



^d 



+ 



+ 



f'^§r + ^o'-i^ + lOr'p^^M] ,p-ii.U 



30 



jj-j 



9 7r2 



a5(/i 



2n2 



1 [ ^^ AW 



30 



1544?) In(^) ^ i /3oln(^)' ^ 1 ^(2)^ 



TT" 



TT 



2 ttAW 



V 



a5(;U ) + ln - 



1\ 5« 



/ 



T^y AW 



'l,8?Bf^.ln(i) 13 fl„Mf^.ln(l) Jl0^i|^ + f %L^^ ___^ ., 



C2 



In i 



re 



5 A(i) 162 TfS A(i) 



30 



TT^ 



+ 



{20 ^tBW loo /3gftA(2) \, /iV 
1 l^¥~4n^ + ^ lO^j^'^i^^j 



30 



TT" 



tt5(/U 



2\5 



16 



/ 



+ 



,PoP.Bi^UnU) , (lo4r + ?^^^^i^)ln(0' l,Pl^'''^4\i\ 



3 7r4 AW 

( 



30 



TT- 



27 7r4 A(i) 



«5(/i')' + 



ttMW ^ V30 7r5 "^6 ttMW 



V 



+ '30 



ttMW V30 
1 30 ze!^ + 30 2^f#^ 7^oln( 



+ - 



^o 



9 TT^ 



i!^i^ I /inE!JoA(2) 



+ — ^ -. ^^ I «5(/i''' 



30 



Tj-o 



i^\ 



+ 



1 30^^|#;^+^5^ ^fMS)_ 

30 TfS + 3 71 



^^(/i^) + 1 



ots{^ 



2n3 



C3 



M^ ^ 160 M^iA(^^ 1^ /^l 



40 ^3^ , 160 ^oftfi^ In 1 ^4 4(2)i„ 1 



30 



+ 



y]-0 



^1 7rM« 



+ s 



/3gfti?'^) , 100 ^?^\,^^ 



1 (40^^^^^^ + fl^^^lln(M 2/5?5(^) 



+ - 



TT" 



3 ttMW 



«5(/" 



2\5 



+ 



+ 



+ 



4/^0 A 5(2) 16'^o^' 'inl - 1 ^ I 4U ^^j + ^ ^^^ I Inl - 



3 Ti^AW 9 TT^AW 30 



TT'^ 



V 



as{^ 



2\4 



/ 



30 7r5 +3 TT^AW 



V 



a5(/i 



2\3 



/ 



1 40 IllfT + 40 %f^ + 40 .» A 20 gglKl)" ! , ,,, 4 ft a.(^^) 



30 



V 



TT"' 



9 7r2 



/ 



3 TT 



^^='30 



50 TfS TfS A(i) 



17 



1 (40 4?? + 120%?^)ln(i 



30 



yj-O 



asifJ, 



2\5 



+ 



1 60 4r + 60-^-afl^ , A?!f^^^ asi.r + S^'^'^^'^^^'''^' 



\ 



30 



n" 



ttMW 



/ 



TT' 



3^(1) 



2\2 



+ 2 



/^0«5(/^^) 



TT^ 



C5 



15 ttMW ^30 ttS ''^^''^ ^ 5 ttMW 



C6 



/5oM(2)«<,(/i2)^ 



TT' 



5^(1) 




Figure 1: DLL A (§) and 6-space (^ results for the transverse momentum distribution 
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Table 1: The first 20 coefficients bm{oo), calculated according to (|T0|). 
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In^i 



are 



61(00) = 62(00) = 64(00) = 0, the coefficients of ( |^ j In ( - ) and 

zero, and for higher orders the biggest power of a In ( - j logarithm is equal to the 



order in 



2tt- 
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Figure 11: The resummed cross section (^), truncated at A^max = 10, compared to 
the results of integration. Top figure: only A^^' coefficient non-zero, middle figure: 
^(1)^^(1) non-zero, bottom figure: A^^\ B^^\ A^'^\ B^"^^ non-zero. The top figure also 
shows a comparison of the effects induced when subsequent coefficients are introduced. 
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Figure 12: Resummation of (0) for the first three orders in ^^(/i^) in the residual 
sum. Here /j^ = Q'^'^qT . Q = Mz = 91.187 GeV, as{Ml) = 0.113. 
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Figure 13: Resummation of ( [21| ) for different subsets of nonzero coefficients, with only 
0(a5(/i^)) terms in the residual sum considered. Here /i^ = Q'^^^qrJ , Q = Mz = 
91.187 GeV, as{Ml) = 0.113. 
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Figure 14: Form factors F'^p^qt), F^i^^qt), with Q = Mz = 91.187GeV, a^lMf) 
0.113. 
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Figure 15: Resummation of the L, NL and NNL 'towers' of logarithms according to (^ 
for n^ = Q2 versus F^i^^qt) from §, with Q = Mz = 91.187 GeV, as{M^) = 0.113. 
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Figure 16: Resummation of the first three (L,NL,NNL) 'towers' according to 0] and the 
first four (L,NL,NNL,NNNL) towers of logarithms (P) for /x^ = q2q^m^^ 91.187 
GeV, as{Ml) = 0.113. 
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Figure 17: Ratio of the results for 3,5,6,7,8 'towers' of logarithms, normalised to the 
4-th tower result, for /x^ = Q^/s^^a^ Q = Mz = 91.187 GeV, as{Ml) = 0.113. 
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